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I give an overview of phenomenological heavy quark results obtained in NRQCD 
on the lattice. In particular I discuss the bottomonium and the b-light hadron 
spectrum. I also review recent results on the decay constants fs and fg 3 . 



1 Introduction 

1.1 Heavy quarks on the lattice 

Lattice QCD provides an approach to calculate the properties of hadronic 
bound states of strongly interacting matter from first principles. However 
when interested in hadronic states involving heavy c and b quarks standard 
lattice methods would lead to large discretisation errors. This is caused by 
the Compton wave length of the heavy quark being non-negligible against the 
lattice spacing. In present simulations the latter is on the order of a few inverse 
GeV. 

However the Compton wave length of the heavy quarks is an irrelevant 
scale for the dynamics of heavy hadrons, see e.g. the lecture note El. One 
possibility to simulate heavy quarks is the expansion of the heavy quark ac- 
tion around its non-relativistic limit, which is known as noii-relativistic QCD 
(NRQCD) SB. Another cure is the heavy clover approaches which is applied 
by several groups presently. Within this talk we will review some of the phe- 
nomenological results obtained in NRQCD. 



1.2 Non-relativistic QCD on the lattice 

In NRQCD the Hamiltonian of the heavy quark is expanded around its non- 
relativistic limit 

D 2 

H = H + 6H , H = . (1) 
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With Hq we denote the leading kinetic term and 5H stands for the relativistic 
and spin dependent corrections. In case of a quarkonium system up to order 
O(MqVq) these corrections read 
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The coefficients c$ have to be determined, such that H in eqn. ([!]) matches the 
Hamiltonian of QCD. So far this has been investigated in lattice perturbation 
theoryB. Since NRQCD is non-renormalisable, the c, will diverge, if the lattice 
spacing a is sent to zero. Therefore the lattice spacing has to be kept finite 
and improvement terms have to be added to the action, such that residual 
discretisation effects become negligible against other sources of error 
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SH disc = c 5 a - c 6 a^-jjj . (3) 



The n is a stability parameter used in equation (Q). Present calculation use 
tadpole improved^ tree-level values for all Cj's. The Hamiltonian (|l|) leads to 
a differential equation being first order in time. Therefore the heavy quark 
propagator Gt can be obtained from an evolution equation 

G t+ i = (1 - a\SH - a\5H Aisc ) (l - a±H )" U+ 

■ (1 - a±H ) n (1 - a\8H - a±5H disc ) G t . (4) 

Thich is numerically quite inexpensive when compared to the matrix inversions 
associated with a Dirac type Hamiltonian. 



2 Heavy quarkonia 

2.1 Bottomonium splittings 

The radial and orbital splittings in the bottomonium (T) and charmonium 
( J /if?) system are « 500 MeV, which approximately equals the average kinetic 
energy of the quarksB. One obtains Vq/c 2 = 0.1 for the T and 0.3 for the J/^f?, 
which justifies the non-relativistic approach. In the case of heavy quarkonia 
the Hamiltonian has to be expanded in powers of vq instead of 1 /Mq 0. 

In figure [I] we show the dependence of spin independent bottomonium 
splittings on the lattice spacing a. This calculation has been performed in the 
quenched approximation, which neglects the vacuum polarisation by quark- 
antiquark pairs. Within the given error bars the lattice result is indeed inde- 
pendent of the lattice spacing. This shows, that within the achieved accuracy 
continuum results can be obtained at finite a. In case_of the ratio of the 
Xb{fP) — Y(1S) energy splitting f^] to the p-meson massa we observe a clear 

a Xb denotes the spin average (xbo + 3xtl + 5xm)/9 
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Figure 1: Scaling of spin independent splittingsE in the quenched approximation. The points 
with prror bars represent the lattice result, the horizontal line the experimental values from 

pdgeI 



missmatch to the experimental result. This is expected to be caused by the 
different running of the strong coupling a s in the quenched theory and the real 
world, where vacuum polarisation is present. This expectation is supported 
by the the ratio of the splittings [T(25) - T( IS) ]/[%;,( LP) - T(15)] being 
closer to the experimental outcome. In this case both quantities probe similar 
momentum scales. 

For spin dependent quantities the equation (Q) is not improved to the same 
level with respect to higher order relativistic and discretisation corrections as 
for spin independent ones. Therefore at this level one observes some scaling 
violation in the bottomonium fine and hyperfine structure El. Also the agree- 
ment of the Xb hue structure with experiment is not as good as for the radial 
and orbital splittings. Improvement needs the inclusion of the higher order 
termsEM and a better determination of the Cj's in the HamiltonianEj. 

2.2 Partly unquenching 

The effect aJLouenching in quarkonium spectroscopy has been investigated by 
two groups ElfiJ. Due to algorithmic reasons, both studies use two flavours of 
dynamical quarks. In figure [2] we compare the outcome for the spin independent 
spectrum to the quenched approximation. Within error bars no significant sea 
quark effect can be shown. 

2.3 Heavy hybrid mesons 

Hybrid mesons denote colour neutral states consisting of gluonic excitations 
as well as a quark-antiquark pair, see e.g. the review by F. Close E3. Of 
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Figure 2: Effect of dynamical fermions on the spin independent T-spectrum. The points 
with error ibars represent the lattice results, the horizontal line the experimental values 
from PDGQ. The a has been determined such that the average of the ^(lP) ~ T(1S) 
and T(2S) — T(IS') splitting fits to experiment. One obtains a - fj ?s 2.4 GeV in all cases. 
The quenched results are taken from the NHQCD collaboration Q. Squares use 2 flavours 
of staggered sea quarks withr«*ri sea = 0.01 Ej, whereas diamonds use Wilson sea quarks 
extrapolated to m aea = m. /3l 16 l. 



special interest are so called exotic states with quantum numbers unavailable 
to gg-mesons, which are hoped to be the first identified experimentally. 

So fa.r .exploratory studies have been undertaken to investigate the bbg 
spectrumliao in NRQCD. These use Hq in eqn. (|4|) only. In this approximation 
one obtains two sets of degenerate states, which we denote by their cubic group 
representation and T 1 ~ + '. States below the Bi3**-threshold are expected 
to be stable. Using the i3}(5732) for this threshold one obtains 11.01 GeV. 
With the a- value from figure [| one gets for the masses 

M{T^~) = 11.09(10) GeV, M(Tf + ) = 11.15(5) GeV . (5) 

The resultEi M(T 1 +_ ) = 10.82(25) GeV is in agreement with the above. With 
this size of an error the question of stability of these states can not be answered. 

3 Heavy light systems 

3.1 Heavy light spectrum 

The physics of heavy light mesons is quite different from the physics of heavy 
quarkonia. For the former in the limit Mq — > oo the heavy quark Q decouples 
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Figure 3: Scaling of B-meson splittings. Please note the different nature of these quantities. 
From the left we give an example for a spin dependent, a flavour dependent and a radial 
splitting. Experimental results apart from the B(2S) are from the PDGQ. The B(2S) is a 
recent result from DELPHI ES 



40 



G 
CP 



a 
a 



400 
, 300 



6.0 

GeV 
5.5 

5.0 



.|.(2^) 



B 



D* 



i(25) 



m) 



B, Bt P-States 



a& Sfc,s^ fib, fi* h 



Figure 4: Speptrum of mesons and baryons containing one heavy 6-quark in file quenched 
approximation EH- Dashed lines indicate experimental results from the PDG □ and dotted 
lines recent results from DELPHI rH, 



from the dynamics and becomes a static colour source. The properties of such 
a state would be determined entirely by the light quarks and the glue. In 
NRQCD this leads to different power counting rules c2l. The corrections in 
eqn. (||) have to be ordered in powers of (AqcdA^q)- 

JSigure || shows results obtained at two different values of the lattice spac- 
ingE2l'E3 in the quenched theory. No sign of large lattice spacing dependence 
can be observed for any of the quantities. For the spin independent results 
one observes reasonable agreement with the experimental result, however the 
results for the i? s -hypernne is in clear disagreement. Here the discussion of 
section 2.1 also applies. The quenched approximation seems to have an ef- 



fect as well, since the heavy clover approach—with different systematic errors, 
delivers a similar reduction of the hyperfineE3. 

Figure |] gives an overview over the hadron spectrum containing one heavy 
6-quark. Apart from the above discussed hyperfine splittings one observes good 
agreement with experimental results, including heavy baryons. 
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Figure 5: Scaling of /b„- Tr^p square@ and octagon^ use higher order corrections in the 
Hamilton than the diamondEJ. The latter has been scaled to match the a- value determined 
from m p used in the other ones. Errors encompass *}«3 statistical uncertainties and those 
arising out of the operator renormalisation and mixingt3. In case of the diamond an estimate 
of the effect of the lower order Hamiltonian is also included. 



3.2 Pseudoscalar decay constant 

Because of the small leptonic branching fraction, the pseudoscalar decay con- 
stant Jb is hard to measure experimentally, however its knowledge is required 
in the determination of e.g. the bag parameter in B-B mixing. 

In order to obtain reliable results the inclusion of operator renormalisation 
and mixing is crucial c3c3'ESc3. In figure |5| we compile a scaling plot for f Bs - 
The plot shows gwd scaling as the lattice spacing is changed. We quote the 
result of referencetJ which is in good agreement with the findings of referenced: 

f b = 147(ll)(i; 2 )(9)(6) MeV, f B . = 175(8)(^ )(11)(7)(±£) MeV . (6) 

These results also jagree with those of other lattice calculations using Dirac 
type Hamiltonians cB. 

4 Conclusion and outlook 

We discussed phenomenological heavy quark results including spectroscopy 
and the decay constants Jb and /b 3 - Future work has to improve on spin 
dependent splittings and quenching effects have to be further addressed. 
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